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The  number  of  test  samples  used  to  characterize  the  fatigue  constants  needed  for 
failure  predictions  for  ceramic  materials  determines  the  confidence  in  these 
predictions.  The  statistical  reproducibility  of  the  dynamic  and  static  fatigue 
experiments  used  to  measure  the  fatigue  constants  was  analyzed  using  both  sta¬ 
tistical  theory  and  a  Monte  Carlo  computer  simulation  technique.  It  was  found 
that  the  statistical  reproducibility  depended  not  only  on  the  number  of  test 
samples  but  also  on  the  other  experimental  test  variables.  It  vies  (continued) 
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Elock  20.  Abstract  (continued) 

shown  that  the  uncertainty  in  the  statistical  reproducibility  can  be  large 
especially  for  sample  size  less  than  about  10C.  Guidelines  for  selecting 
the  optimum  sample  size  for  a  given  dynamic  or  static  fatigue  experiment  are 
given.  It  is  recommended  that  before  meaningful  conclusions  can  be  drawn 
regarding  the  effect  of  a  test  variable  on  fatioue,  the  statistical  reproduc¬ 
ibility  of  tne  experiment  be  determined. 
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FOREWORD 


This  report  describes  the  results  of  an  experimental  program 
oriented  toward  a  better  understanding  of  lifetime  predictions  for 
optical  glass  fibers.  Some  of  the  progress  made  toward  this  goal 
is  summarized  in  the  attached  technical  paper  comprising  this 
report. 
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ABSTRACT 

The  number  of  test  samples  used  to  characterize  the  fatigue  constant: 
needed  for  failure  predictions  for  ceramic  materials  determines  the  con¬ 
fidence  in  these  predictions.  The  statistical  reproducibility  of  tne 
dynamic  and  static  fatigue  experiments  used  to  measure  the  fatigue 
constants  was  analyzed  using  both  statistical  theory  and  a  Monte  Carlo 
computer  simulation  technique.  It  was  found  that  the  statistical  reproduci¬ 
bility  depended  not  only  on  the  number  of  test  samples  but  also  on  the  other 
experimental  test  variables.  It  was  shown  that  the  uncertainty  in  the  statistical 
reproduc ibil ity  can  be  large  especially  for  sample  size  less  than  about 
100.  Guidelines  for  selecting  the  optimum  sample  size  for  a  given  dynamic 
or  static  fatigue  experiment  are  given.  It  is  recommended  that  before 
meaningful  conclusions  can  be  drawn  regarding  the  effect  of  a  test 
variable  on  fatigue,  the  statistical  reproducibility  of  the  experiment 


be  determined. 


1.0  INTRODUCTION 


Methods  g*  dealing  with  design  problems  involving  fatigue  of  ceramic 
materia  . s  neve  been  developed  over  the  past  10  years  through  tne  application' 
of  fracture  mechanics  principles.  Since  these  principles  can  be  used  to 
characterize  both  J»e  conditions  for  subcritical  crack  growth  anc  the 
conditions  for  crack  instability,  they  can  be  used  for  purposes  of  design, 
to  estimate  the  allowable  stress^or  the  expected  lifetime,  or  trie  proc* 
stress  necessary  to  assure  a  minimum  lifetime.  This  is  accomplished  by 
estimating  the  initial  crack  size  in  a  ceramic  component  end  tne  time 
required  for  this  initial  crack  to  grow  to  a  critical  size  for  spontaneous 
fracture.  For  example,  it  has  been  derived  by  assuming  a  simple  power  law 
relationship  between  subcritical  crack  velocity  and  stress  intensity  tr.ct 
the  failure  time  (t,)  under  a  constant  applied  stress  (c  )  is:"'  * ^ 

•  d 


r  B 


2  N-2 

where  B  =  2/(AY  (N-2)  K.^.  ),  A,  N  =  material/environment  constants,  Y  = 

geometric  constant  (about  1.2  for  surface  flaws),  =  critical  stress 
intensity  facto",  and  *  fracture  strength  in  an  inert  environment.  B  and 
N  in  Eq.  (1)  are  fatigue  constants  that  for  a  given  material /environment  system 
characterize  subcritical  crack  growth.  The  inert  strength  in  Eq.  (1)  character¬ 
izes  the  initial  flaw  size.  If  p'oof  testina  is  used  to  truncate  the  flaw 
distribution,  then  the  minimum  inert  strength  after  proof  testing  is  equal 
to  the  maximum  proof  stress  [rr),  hence,  tne  correspondingly  minimum  failure  time 


Kl  -  3- 


Froiii  sq.  .1  and  {,2}  it  i  s  seen  that  failure  t.-red  iclimit  *r«-  neienaer 
oil  the  Tutique  | ki rameters  N  arid  ii.  These  pai  amutfi  art.-  v  mi-, lair.  f t»r  a 
given  materiu  ! /env  t  roriment  system  and  <_an  be  experimentally  deter::  \  ned 

3 

directly  using  fracture  mechanics  techniques,  or  can  be  indirectly  measured 

2 

using  static  fatigue  or  dynamic  fatigue  experiments.  unfortunately,  failure 

predictions  are  extremely  sensitive  to  the  experimental  uncertainty  it;  the 

fatigue  parameters.  Statistical  techniques  for  estimating  tins  uncertainty 

4  5 

in  failure  predictions  have  been  developed,  '  however,  the  statistical 
reproducibility  of  the  experimental  techniques  used  to  evaluate  the  fatigue- 
parameters  !i  and  5  has  not  been  previously  determined.  Statistical  repro¬ 
ducibility  is  due  to  random  sampling  errors  that  are  inherent  in  every  ex¬ 
periment.  In  the  random  selection  of  a  finite  number  of  samples  for  testing, 
one  would  expect  to  see  some  statistical  variability  in  the  measured  properties 
that  would  be  dependent  on  the  number  of  samples  selected.  It  further  would 
be  expected  that  this  variability  in  the  estimation  of  the  fatigue  parameters 
would  increase  as  sample  size  decreases;  however,  this  statistical  reproduc i bi 1 i 
has  never  been  quantified  although  some  Monte  Carlo  results  on  the  re¬ 
producibility  of  N  as  determined  by  dynamic  fatigue  tests  have  been  previously 
reported  by  the  present  authors.^  Before  meaningful  conclusions  can  be  drawn 
from  the  results  of  fatigue  experiments,  the  statistical  reproduc i b i 1 i ty  of 
these  experiments  must  be  known. 

The  purpose  of  this  paper  is  to  quantify  the  statistical  reproducibility 
of  the  dynamic  and  static  fatigue  experiments  for  measuring  the  fatigue 
constants  N  and  5.  Statistical  reproducibility  is  analyzed  using  both 


statistical  theory  and  a  Monte  Carlo  computer  simulation  technique.  Since 


the  statistical  theory  contains  a  number  of  critical  assumptions,  it  is 
important  to  independently  validate  the  statistical  theory  approach  with 
tne  Monte  Carlo  techr.iaue.  Emphasis  is  placed  on  tne  dynamic  and  static 
fatigue  tecnniques  because  they  are  increasingly  being  used  to  measure 
the  fatigue  behavior  of  ceramics.  This  is  because  these  test  techniques 
can  util  ice  samples  containing  flaws  representati ve  of  those  on  actual 
components.  It  is  believed  that  the  results  of  this  study  will  leaG  to  a 
better  understanding  of  the  static  and  dynamic  fatigue  test  techniques  and 
their  statistical  variability.  From  this  information  guidelines  for  sample 
size  reauirements  in  terms  of  optimum  statistical  reproducibility  can  be 
developed. 


I 


2.0  ANALYSIS  OF  STATISTICAL  RE PF.ODUC I E I L I TY 


2.1  Dynamic  Fatigue 

Dynamic  fatigue  data  is  generated  by  measuring  the  fracture  strength  of 

a  number  of  samoles  at  several  constant  stressing  rates.  The  fatigue  constants 

N  and  E  can  be  determined  from  dynamic  fatigue  data  tnrougn  using  one  ox  four 

analyses:  median,  homologous  stress,  iterative  bivariant,  or  iterative  tri- 
2  7 

variant.  ’  Since  all  of  these  techniques  analyze  the  same  set  of  dynamic 
fatigue  data  in  determining  N  and  E  through  a  linear  regression  analysis,  all 
are  expected  to  result  in  essentially  the  same  statistical  reproducibility 
for  N  and  B.  Thus,  the  median  analysis  technique  was  chosen  for  this  study 
because  of  its  simplicity  and  wide  usage. 

2  7 

With  the  median  analysis  the  dynamic  fatigue  data  are  fitted  so:  ’ 


r. 


(3) 


where  S  =  median  fatigue  fracture  strength,  c  =  stressing  rate,  and  a,,  a?  - 

°  2 
linear  regression  constants.  The  fatigue  constants  are  then  determined  from:  ’ 


'  a, 

„  /)  /  ^  A 

>'r.  -  N  i  I  )<X.  -  ■  fJ*i)  -  f  fj .  i'. 

"  ,'j  i  -y-**  1  "  i  7  3  ;  £  ■ 

Wl  ,  l 

where  =  median  inert  strength. 

Dynamic  fatigue  strength  data  of  an  "ideal"  material  was  simulated 

8  0 

on  a  computer  using  a  Monte  Carlo  technique.  It  was  assumed  that  the 


(4a) 

(4b) 


fatigue  constants  N  and  B  of  this  ideal  material  are  given  and  that  the 
inert  strength  distribution  is  given  by  a  two  parameter  Weibull  distribution 


whose  slope  and  scale  parameters,  riband  $0 ■  respectively,  are  known.  With 
the  Monte  Carlo  technique  a  given  number  of  samples  at  a  specific  stress¬ 
ing  rate  were  chosen  by  randomly  selecting  their  failure  probati’ity 
from  a  umforn  distribution  between  C  and  1.  The  corresponding  *'atigue 

fracture  strengths  were  then  calculated  according  to  the  fracture  mechanics 
2 

relationship: 

l  C;1  r  ^  ,77r  +  *  J-  <r  J 

aj*  t 

A  separate  set  of  inert  strength  samples  were  chosen  similarly  to  the  fatigue 

samples,  namely  failure  probabilities  were  randomly  chosen  and  tne  corresponding 

1  0 

inert  strength  were  calculated  from  the  two  parameter  WeiPull  distribution: 


Once  a  set  of  inert  strength  and  fatigue  strengths  at  several  different 
stressing  rates  were  randomly  chosen,  N  and  B  were  determined  from  this  data 
using  Eqs.  (3)  and  (4).  It  should  be  noted  that  the  median  strength  for  a 

given  set  of  strength  values  was  determined  by  ranking  and  fitting  the  strengths 
to  a  Weibull  distribution  by  linear  least  square  method.  The  meaiar.  strength 

was  then  calculated  as  the  value  at  F=  0.50.  Alternatively,  tne  median 
strength  could  have  been  determined  by  choosing  the  actual  median  strength 
value;  however,  this  results  in  greater  variability  and  thus  was  not  used.  By 
iterating  this  procedure  100  times,  distributions  for  N  and  B  were  generated 
which  reoresent  the  statistical  reproducibility  of  N  and  B  as  determined 
from  the  dynamic  fatigue  test.  With  the  100  values  generated  by  the  Monte 
Carlo  technique,  the  average  values  with  their  corresponding  variances  and  the 
covariance  between  N  and  B  were  calculated  from  the  usual  statistical  formulas:'1'1 


/J  ' 


5  /j 


J  O  ~  n 


1/^6)  '  “  £  ^-3,  -  -^S  ) 


1  'n''  J'» 


where  n  =  100  for  these  computer  calculations.  For  a  given  material/environment  system, 

i.e.  for  a  given  set  of  m,  S  ,  K,  and  B,  the  important  test  variables  studied  were 

o 

the  number  of  samples  per  stressing  rate,  the  stressing  rate  range,  and  the 
number  of  stressing  rates  that  were  used  in  a  given  test.  For  the  case  of 
mul ti -stressing  rates,  the  stressing  rates  were  evenly  spaced  from  maximum 
to  minimum.  Figure  1  gives  a  schematic  flow  diagram  of  the  Monte  Carlo 
computer  simulation  technique  for  determining  the  fatigue  constants  N  and  B 
by  the  dynamic  fatigue  test.  Note  that  the  same  number  of  inert  samples  were 
chosen  as  the  number  of  samples  per  stressing  rate. 


i  he  variances  and  covariance  of  the  inert  sc>"enctr:  parameters  we^e  a' 
derived  as  a  function  of  sample  sice  from  statistical  theory  ov  makir.r 


certain  assumptions  (see  Appendix,  Eos.  (Al£) ,  < A26) ,  anc  : 


H  C  i  1  .  i  Ti  6  S  L 


of  this  analysis  are  as  follows: 
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where  J  =  number  of  samples  and  f  =  gamma  function.  Likewise,  the  variances 


the  fatigue  constants  K  and  E  and  their  covariance  were  deriv'ec  as 


s  nown  i n  t 


ppendix ,  Eos.  (A41),  (A42),  and  (MS),  to  be: 
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where  Jq  =  total  number  of  samples  used  in_  dynamic  fatigue  experiment 
number  of  stressing  rates,  f\(S C  .  =  ~  Z  ( <T.  -  *-*■.  <r  )  /  d 
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Since  these  derivations  made  a  number  cf  critical  assumptions,  it  is  o* 
importance  to  compare  these  variances  witn  tnose  genera  tec  by  tne  Home  Carle 
teennique  to  aeterrrine  tne  validity  of  cne  assumptions. 

2.2  Scatic  Fatigue 

Static  ~atigue  tests  entail  the  repeatec  measurement  o*  *a''lure  time  at 

several  constant  app’iiec  stresses.  Static  fatigue  oata  can  oe  analyzed 

similar  to  dynamic  fatigue  data  by  the  median,  homoloaous  stress,  iterative 

2  7 

bivariant,  and  iterative  trivariant  analyses.  ’  As  in  the  case  cf  dynamic 

fatigue  the  median  analysis  was  chosen  for  studying  the  statistical  repro- 

2  7 

ducibilitv.  With  tne  median  ancivsis  the  data  is  fitted  to: 


t  _ 


(1C) 


where  t,  =  median  failure  time  and  a.,  a,  =  linear  recression  constants, 
f  j  4 

2 

The  fatiaue  constants  N  and  B  are  determined  from:  ’ 


(He) 


(lib) 


The  Monte  Carlo  analysis  of  the  statistical  reproducibility  of  A  and  B 

as  determined  from  static  fatigue  data  is  similar  to  that  used  for  dynamic 

fatigue  data.  First,  a  given  number  of  failure  times  at  a  specific  applied 

stress  are  selected  by  randomly  choosing  the  failure  probability  and  then 

2 

calculating  the  corresponding  failure  time  from: 
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r, 

A  r- 
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CsJ 


Once  a  set  of  failure  times  at  several  different  appliec  stresses  and  a  set 
cf  inert  strengths  are  randomly  chosen,  N  and  E  are  determined  fro;;.  Ecs.  ’IE'1 


1 


and  (11;.  By  iterating  tnis  procedure  100  times,  a  distribution  cf  U  and  E 
values  are  generated  from,  which  the  statistical  reproauci bil i ty  o*  these 
parameters  can  be  determined.  Similar  to  the  case  of  dynamic  *atigue,  the 
variables  studied  for  a  given  material/envi ronment  system  were  the  number  cf 
samples  Der  applied  stress,  the  applied  stress  range,  and  the  nui-oer  cf 
applied  stresses  in  a  given  test. 

The  variances  in  the  fatigue  constants  N  and  B  and  their  covariance  co„lc 
also  be  derived  from  statistical  theory  as  shown  in  the  Aooencix  ,  Egs.  (A53  ) , 


( A 5 3 ) ,  and  (A54)  to  be: 
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where  now  J  =  number  of  samples  used  in  static  fatioue  experiment,  J„  =  number 
o 

of  applied  stresses,  =  number  of  samples  tested  at  each_applied  stress, 
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Again  it  was  of  interest  to  compare  these  variances  to  these  determined  by  tne 
Monte  Carlo  technique  to  determine  the  validity  of  the  assumptions  made  in  ceri' 
inc  the  above  equations. 


3.C  RESULTS  AND  DISCUSSION 


Both  the  Weibull  parameters  (m  anc  S  )  and  the  fatigue  constants  (I*  anc 

depend  stronclv  on  sample  size.  Figures  2  and  3  snow  the  dependency  for  n  an 

S  .  respectively  on  sample  size  where  the  coefficient  cf  variation  is  definec 
c 

as  the  standard  deviation  divided  by  the  mean  value  of  tne  parameter  in 
question.  From  tnese  figures  it  can  be  seen  that  there  is  good  agreement 
between  the  Monte  Carlo  technique  for  estimating  the  statistical  reprodu¬ 
cibility  and  that  derived  from  statistical  theory.  Tne  figures  further 
show  that  the  statistical  variability  in  the  Weibull  slope  parameter  m  is 
a  function  of  only  sample  size  while,  for  the  Weibull  scale  parameter  S  , 
the  variability  is  a  function  of  both  sample  size  and  the  Weibull  slope. 

Since  for  small  sample  sizes  the  statistical  uncertainty  in  the  Weibull 
parameters  can  be  large,  especially  for  m,  important  judgements  and 
significant  analyses  of  strength  should  not  be  based  on  smell  sample  sizes. 
From  Fig.  2  ana  3  it  can  be  seen  that  sample  sizes  of  at  least  30  should  be 
used  for  all  but  the  most  preliminary  investigations,  although  for  small  m's 


tc  get  acceptable  levels  of  So  may  require  as  many  as  100  samples.  Similar 
results  and  conclusions  were  reached  in  earlier  Monte  Carlo  simulation 


studies . 


12,  13 


Figures  4  and  5  show  that  the  statistical  reproduci bi 1 i ty  cf  the  fatigue 
parameter  N  and  B  as  determined  by  the  dynamic  fatigue  experiment  is  strongly 
dependent  on  both  sample  size  and  the  fatigue  resistance  of  the  materia1 
(large  N  values  generally  represent  materials  with  a  greater  fatigue 
resistance).  For  sample  sizes  of  less  than  100,  the  statistical  uncertainty 


in  N  ana  B  can  be  very  large,  especially  for  the  more  fatigue  resistant 
material.  Again  there  is  good  agreement  between  the  Monte  Carlo  technioue  fo 


theory;  thus,  giving  evidence  of  tne  validity  of  the  assumptions  made  in 
deriving  the  equations  for  statistical  variability  (see  Appendix). 

Tne  statistical  variability  of  the  fatigue  parameters  !i  and  B  as  determine 
by  tne  dynamic  fatigue  experiment  was  also  dependent  on  the  number  cr  stressing 
rates  chosen,  stressing  rate  range,  and  the  Weibull  slope  rr. .  For  examrle. 

Fig.  6  shows  that  for  the  same  range  of  stressing  rates  (maximum  to  minimum; 
and  the  same  total  number  of  samples,  uncertainty  in  statistical  reproduci bi 1 i t 
increased  as  the  number  of  stressing  rates  used  in  determining  !-.  is  increased 
from  2  to  7.  The  best  statistical  reproducibi 1 ity  occurs  for  the  case  where 
N  is  determined  from  strength  measurements  at  two  stressing  rates  corresDondi nc 
to  the  maximum  and  minimum.  Figure  7  shows  the  statistical  variability  of 
as  a  function  of  m,  keeping  the  other  parameters  constant,  and  illustrates 
that  low  m  values,  corresponding  to  a  greater  variability  in  strength,  re¬ 
sult  in  a  larger  uncertainty  in  N.  Figure  8  shows  that  the  statistical 
variability  of  K  is  quite  sensitive  to  the  stressing  rate  range,  with  the 
variability  increasing  with  a  decreasing  range  of  stressing  rates.  Finally, 
it  should  be  noted  that  the  Monte  Carlo  results  were  left  out  of  Figs.  6,  7, 
and  8  for  clarity;  however,  these  results  agreed  quite  well  with  those  shown  in 
the  figures  based  on  statistical  theory. 

The  distributions  of  the  fatigue  parameters  N  and  B  as  generated  by  the 
Monte  Carlo  computer  simulation  technique  for  the  dynamic  fatigue  experiment 
could  be  approximated  by  a  normal  distribution.  Finure  ?  shows  that  a  normal 
distribution  well  represents  the  histogram  for  the  fatigue  parameter  N. 

Histoarams  for  the  Weibull  parameters  m  and  S  could  also  be  approximated  bv 

o 

normal  distributions.  This  is  important  because  the  confidence  limits  of 


reproducibility  for  a  given  parameter  can  be  estimated  by  simply  multiplying 


r.e  stancaro  deviation  bv  me  aocrocriace  'actor  -or  nc-ma,  oiscr-r;. 


’us ,  one  35”  reprocucibi lit/  1  ini ts  would  be  -1 . 360  stancaro  oev a tiers ’  arc 


-eoreser.t  ore  interval  .-mere  aoout  rob  or  one  va.ues  -or  a  gi  /en  cara-reae-  .  *a. 
“r.e  3Cb  reorccucibil lay  1  imias  .-/cult  be  -1 .  5*2-5  stancaro  oevaocors  1  arc  so  or. 

Results  on  ore  staoi so: cal  reorpduci bi 1 i ay  b"  ore  -atigue  oaraoeoars 
.‘i  arc  3  as  determined  by  the  static  fatigue  experiment  were  similar  oo 
ones  a  shewn  for  the  dynamic  fatigue  experiment.  In  particular,  t.re  statistical 
variability  of  '<  anc  3  as  determined  oy  ore  static  -atigue  exoerimeno  decreased 
with  increasing  sample  size,  decreased  for  material s  with  lower  '!  values, 
decreased  as  ore  aooliec  stress  -ance  increased,  and  decreased  as  ore 
we i bull  sicca  oa ■•ameter  m  increasac.  "or  the  same  applied  stress  range  anc 
nur.cer  o*  total  samples,  one  best  reorocuci bi  1  ’  ty  occurs  for  the  case  where 
.1  ana  3  are  ceterai ned  ~r cm  time-to-fai lure  measurements  at  two  acoliec 


;sses  corresponding  to  the  maximum  anc  minimum.  For  t re  case  o-  mul  v-asci  'sc 


stresses,  the  statistical  variability  depended  mainly  on  the  total  number  of 
samDles,  not  on  the  number  of  samples  per  applied  stress,  and  only  somewnat 
on  the  number  of  apolied  stresses  used  to  determine  N  and  3.  As  before, 


there  was  gocc  agreement  between  t.ne  Monte  Carlo  analysis  and  t.ne  stato  stoca  i \-j 


derived  variability  equations  ,see  ,-opencix : .  A, so,  t.ne  distributions  generated 


by  the  Monte  Carlo  analysis  could  be  approximated  by  normal  distributions. 
Finally,  io  should  be  noted  that  the  magnitude  of  t.ne  statistical  reprocucibi 
of  'I  and  3  determined  by  t.ne  stat’d  fatigue  exoer-'ment  is  s’.mi"!ar*  oo  t.nat 
:e terr.ir.ee  *rom  t.re  tvram’c  -atioue  excer-nent  for  avoids"  -areas  of  c/namic 


-V 


1 


h 
\ ' 


j 
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Finally,  it  should  be  noted  that  there  was  a  small  bias  present  in  the 
Monte  Carlo  generated  values  for  m,  S  ,  N,  and  B,  especially  for  small 
sample  sizes.  It  is  well  known''  that  when  random  variables  are  combined  in  a 
non-linear  fashion,  the  resulting  quantity  is  generally  subject  to  biases,  i.e. 
the  combined  effect  of  the  random  fluctuations  of  the  individual  variables 
will  cause  the  derived  quantity  to  be  systematical ly  larger  or  smaller  than 
it  would  have  been  in  the  total  absence  of  such  fluctuations.  Biases  in  m 
and  Sq  have  been  previously  discussed.  ’ ' J ’  3  These  systematic  biases  in  the 
derived  parameters  were  quite  small  in  comparison  with  the  statistical 
reproducibility  of  the  parameters;  hence,  they  were  neglected  in  the 
statistical  reproduci bi 1 i ty  of  a  given  parameter. 


.  0  . 


practical  :cr-!s::ESAT:o:-is 

Cnee  t.ne  s ta "1  szica'  '•eorccuc'bility  o~  tr.e  ~ a ~  1  e: u e  car abaters  'i  ere  £ 

■  s  known  as  e  "unction  g-  sate’e  site  ar.c  t re  other  test  /ar* aoi es ,  the 
corresccr.ci  nc  statistical  uncertainty  in  "ai’ure  creci eti ers  car  t r.er.  ce 
cecer-iinec  as  ciscussec  zre-r.  ;esi  j.  ”  ’  *  Is  -3  -rzr, i  t.nis  knew! ecce  era:  tne 
numcer  of  sane'es  recu’re-c  "or  2  given  oynamic  or  static  "argue  test  can 
be  determiner  ~'ne  exact  cro’ce  o~  same'e  size  *cr  a  riven  test  wiV.  cece°d 
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in  the  dynamic  fatigue  test.  An  uncertainty  of  +  10*  requires  about  140 
samples  to  be  used  in  the  dynamic  fatigue  test  for  the  stressing  rate  range  of 
0.013  to  50  MPa/s.  However,  it  can  be  seen  from  the  Fig.  12  that  if  the  stress¬ 
ing  rate  range  is  increased  to  0.0005  to  50  MPa/s,  then  the  number  of  samples 
required  woulo  decrease  to  about  75  On  the  other  hand,  if  the  stressing 
rate  range  is  decreased  to  0.05  to  5.0  MPa/s,  the  number  of  samples  required 
increases  to  about  300.  This  figure  clearly  shows  that  sample  size  require¬ 
ments  are  not  necessarily  small  and  depend  on  the  specific  degree  of  reproduci¬ 
bility  that  is  acceptable. 

The  analysis  of  statistical  reDroduci bi 1 i ty  can  also  be  useful  in  de¬ 
termining  whether  a  given  variable  effects  the  fatigue  benavior  of  a 
material.  For  example,  the  fatigue  behavior  of  soda-lime  glass  has  been 
measured  as  a  function  of  test  environment  (6N  NaOH,  distilled  water,  and 
6N  HC 1 )  using  the  dynamic  fatigue  test  with  stressing  rate  range  of  « 6  to 
0.17  MPa/s  and  a  total  of  120  samples  for  each  dynamic  fatigue  test.16 
The  effect  of  test  environment  can  be  seen  by  calculating  the  allowable  stress 
for  samples  that  have  been  proof  tested  up  to  ICO  MPa  and  that  must  survive 
a  minimum  of  one  year  in  service.  The  result  of  these  calculations  are  given 
in  Table  I  with  the  indicated  +  one  standard  deviation  limits  due  to  the 
statistical  uncertainty  in  the  fatigue  constants  N  and  B.  It  is  seen  that 
the  three  predicted  allowable  stresses  are  well  outside  of  a  standard  deviation 
from  each  other  and,  thus,  it  is  likely  that  the  test  environments  do  have  a 
significant  effect  on  the  fatigue  behavior  of  soda-lime  glass.  While  this 
conclusion  seems  quite  straightforward,  it  must  be  remembered  that  before  any 
meaningful  conclusions  can  be  drawn  regarding  the  effect  of  a  particular 
variable  on  the  fatigue  behavior  of  a  material,  it  must  be  demonstrated  tnat 
the  effect  is  larger  than  the  statistical  reproducibil ity  of  the  experiment. 
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Table  1.  Predicted  Allowable  Applied  Stress  for  a  Minimum  Lifetime  of 
One-Year  for  Soda  Lime  Glass  Samples  Proof  Tested  at  100  MPa 


Environment  _ N _ 1  n  3  (MPa^ •  s  )  va 


6N 

NaOH 

13. 

,5 

(+  2.31) 

0. 

.163 

(± 

1.64) 

25. 

.35 

(l 

3. 

.36) 

Dis 

tilled  H^O 

13. 

.0 

(+  1.11) 

2 , 

.585 

(t 

0.37) 

1  5 . 

.32 

(± 

1  . 

.53) 

6N 

HC1 

25. 

,1 

(+  4.35) 

-16. 

.053 

(  + 

4.67) 

13, 

.36 

(  + 

2 , 

.15) 

+  Number  in  parenthesis  is  +  one  standard  deviation  corresponding  to  statistical 
reproduci bi 1 i ty  of  dynamic  fatigue  experiment. 
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Figure  2 

Figure  3 

Figure  4. 

-igure  5 

Figure  5 

Ficure 


Flow  aiagran  for  Monte  Carlo  tecnniaue  for  evaluating  the  statistical 
reoroducibi  1  i ty  of  "<  ana  3  as  determined  from  the  dynamic  fatigue 
experiment . 


"he  coefficient  of  variation  (C .'■/.}  of  Weibuil  sloce  parameter  m 
as  a  function  of  sample  size  for  various  values  of  m.  The  rente 
Carlo  results  are  given  by  the  data  Qoints  and  are  compared  to 
predictions  from  statistical  theory,  Eq. (3a). 

The  coefficient  of  variation  of  Weibuil  scale  parameter  S  for 
various  values  of  m  as  a  function  of  sample  size.  The  Moflte  Carlo 
results  are  given  by  the  data  points  and  are  compared  to  predictions 
from  statistical  theory,  Ea.l3b). 

The  coefficient  of  variation  of  fatigue  parameter  1!  as  a  function 
of  total  number  of  samples  used  in  dynamic  fatigue  experiment  with 
3  stressing  rates  (0.313  to  50  MPa/s).  The  Monte  Carlo  results  are 
given  by  the  data  points  and  are  compared  to  predictions  from  statisti¬ 
cal  theory,  Ea.(3a'.  Th^  Weibuil  and  fatigue  constants  chosen  are 
typical  of  soda-lime  glass-  (m.=  15.0,  S  .=  133.0  MPa,  11  =  13.4,  and 

?  '  C  1  n  j 

3  =  0.13  MPa"-s)  and  vitrified  Grinding  wheel  material1 
'm^=  15,  S  .*  53  MPa,  M  =  <13.2  and  3  =  1.1  x  10°  MPa"-S/. 

The  coefficient  of  variation  of  fatigue  parameter  5  as  a  -'unction  of 
total  number  of  samples  used  in  dynamic  fatigue  experiment  with  3 
stressing  rates  (0.013  to  50  MPa/s).  The  Monte  Carlo  results  are 
given  oy  trie  sata  points  and  compared  to  predictions  -rom  statistical 
theory,  Eg.(ab).  "he  Weibuil  and  fatigue  constants  chosen  are  typical 

of  soda-lime  giass^  (m.=  1 3 .  0, 3 _ . =  133.0  MPa,  N  =  13.4,  and  3  = 

n  1  31  1  i 

0.13  MPa"-5/  and  vitrified  grinding  wheel  material'  ;m,=  15, 

3, .  =  53  '-'Pa,  .'I  =  43.2  and  3  =  1  . 1  x  ICf3  MPa^-s). 

O  1 

The  coefficient  of  variation  of  fatigue  parameter  '!  as  a  -‘unction  of 

total  number  of  samples  used  in  dynamic  fatigue  experiment  with 

stressing  rate  range  of  0.313  to  50  MPa/s.  The  number  of  stressing 
rates  used  in  determining  4  is  as  indicated  and  the  appropriate  constants 

are  11  =  13.4,  3  =  3.13  MPa’s,  n.=  3.2,  and  S  .=  133.0  MPa. 

1  o  i 

"he  coefficient  of  variation  of  fatigue  parameter  U  as  a  -’unction  of 

total  number  of  samples  used  in  dynamic  -aticue  exoer'ment  -'or  3 

stressing  rates  (3.313  to  30  MPa/.s)  ana  as  a  ‘‘unction  of  Weibuil 
3 1  ooe  parameter  it .,  while  tne  ptner  constants  a^e  'I  =  13.- 

3  =  0.13  MPa~>s,  and  5Q--=  130.0  MPa) . 


Figure  3.  The  coefficient  of  variation  of  fatigue  parameter  "  as  a  *uncti;n  of 

total  number  of  samples  used  in  dynamic  fatigue  experiment  -'or  3 

stressing  rates  with  two  ranges :  0.01 3  to  50  ana  0.05  to  5.0  3Ma/s. 

The  appropriate  constants  are  '1  =  13.4,  3  =  0.13  MPa'-s,  5.2, 

ana  S01-=  138.0  MPa. 

Figure  9.  Histogram  for  fatigue  parameter  'I  as  determined  Monte  Carlo  simulated 
dynamic  fatigue  experiment  for  3  stressing  rates  ( C .013  to  50  VPa/3; 
using  a  total  of  150  samples.  The  superimposed  smooth  curve  is  the 
probability  density  function  of  the  normal  distribution  with  t re  same 
mean  and  variance  as  the  Monte  Carlo  results. 

Figure  10.  The  coefficient  of  fatigue  parameter  .'!  as  a  function  of  total  number 

of  samples  used  in  dynamic  fatigue  experiment  with  3  stressing  •'ates 

(0.013  to  50  MPa/s)  and  static  fatigue  exoeriment  with  3  applied  stresses 

7 

(40  to  53.6  MPa).  Experiments  are  typical  for  soda-lime  glass" 

(11  =  13.4,  B  =  0.13  MPa2‘S,  1^.=  3.2,  and  S  ^  133. 0  MPa). 

Figure  11.  Design  diagram  for  minimum  lifetime  prediction  after  proof  testing  rcr 

soda  lime  glass  (M  =  13.4,  3  =  0.18  MPa2*s,  m, =  3.2,  and  S  .=  133. C  "Pa) 

*  0  i 

The  90‘3  confidence  limits  are  shown  for  3  sample  sizes  (45  ,  20,  and  130} 
used  in  dynamic  fatigue  experiment  with  3  stressing  rates  (0.313  to 
50  MPa/s). 

cigure  12.  90"  confidence  limits  for  the  proof  stress  to  assure  a  minimum  lifetime 
of  10  years  at  a  constant  applied  stress  of  7  MPa  for  soda-lime  glass 

(11  =  13.4,  3  =  0.18  MPa2-s,  m;=  3.2,  S  .  =  133.0  MPa)  as  a  -'unction 
of  sample  size  used  in  dynamic  fatigue  eSberiment  with  3  stressing  rates 
in  the  range  of  0.05  to  5.0,  0.013  to  5C,  and  0.005  to  50  MPa/s. 
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A. 1.  Linear  Regression 

Linear  regression  is  generally  usee  to  analyze  one  test  data  ir.  rietfcrr-'n- 

inc  one  fatigue  parameters  N  and  B.  Trie  variances  of  these  fatigue  oa-amete**' 

can  chen  be  estimated  ov  using  cne  Ian  cf  propagation  of  errors .  commonly 

11  i ' 

known  as  the  chain  ruie.  . he  chain  rule  is  a  linear  approximation  and 

becomes  inaccurate  when  the  coefficient  of  variation  (ratio  of  the  standard 
deviation  to  tne  mean'  of  tne  data  is  large.  Manael  ^  !  gives  the  rule  of  thumb 
that  tne  coefficient  of  variation  should  not  exceed  10L. 

In  linear  regression  J  data  values  are  least  squares  fitted  to  a  straight 
line  wncse  equation  is  given  bv: 

o  -  dL  ?L  +  k  :  j_ 


where  y  is  the  dependent  variable,  x  is  the  indeDendent  variable,  a  is  the  slope 
and  b  is  the  intercept.  Assuming  that  all  the  inaccuracies  occur  in  the 
measured  values  of  y  and  that  tne  parent  standarc  deviation  of  the  data,  r, 
is  independent  of  x,  the  variances  and  covariance  of  the  parameters  a  and  b 
can  be  expressed  as : 1 1 
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In  Eqs.  (A52)-(A54)  it  is  assumed  that  the  covariances  between  a3,  a^  and 

A 

In  S  are  zero  and  that  the  number  of  inert  strength  samples  is  the  same  as  the 
number  of  samples  tested  in  fatigue  at  each  applied  stress. 
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